Statistics of the Coulomb blockade peak spacings of a silicon quantum dot 
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We present an experimental study of the fluctuations of 
Coulomb blockade peak positions of a quantum dot. The dot 
is defined by patterning the two-dimensional electron gas of 
a silicon MOSFET structure using stacked gates. This per- 
mits variation of the number of electrons on the quantum dot 
without significant shape distortion. The ratio of charging 
energy to single particle energy is considerably larger than 
in comparable GaAs/AlGaAs quantum dots. The statistical 
distribution of the conductance peak spacings in the Coulomb 
blockade regime was found to be unimodal and does not fol- 
low the Wigner surmise. The fluctuations of the spacings 
are much larger than the typical single particle level spacing 
and thus clearly contradict the expectation of random matrix 
theory. 

PACS numbers: 73.23.Hk,05.45.+b,73.20.Dx 



The spectral properties of many quantum mechanical 
systems whose classical behavior is known to be chaotic 
are remarkably well described by the theory of random 
matrices (RMT) |Q. This has been experimentally con- 
firmed, for example, in measurements of slow neutron 
resonances of nuclei Q and in microwave reflection spec- 
tra of billiard shaped cavities Q . Electron transport ex- 
periments performed on semiconductor quantum dots in 
the Coulomb blockade (CB) regime Q provide a further 
possibility to check RMT predictions. The classical mo- 
tion of electrons in these structures can be assumed to be 
chaotic due to an irregular potential landscape produced 
by impurities, an asymmetric confinement potential ||, 
and/or electron-electron interactions ||. The transport 
properties of quantum dots are inherently related to their 
energy spectra and electronic wavcfunctions and thus the 
connection with RMT is readily made ||[7) . 

Indeed, experiments on the distribution of conductance 
peak heights of quantum dots in the Coulomb blockade 
regime have shown good agreement with the predictions 
of RMT HI). On the other hand, the distribution of 
the CB peak spacings have been found to deviate from 
the expectations of RMT 10-12J. The results suggest 
that the peak spacings are not distributed according to 
the famous Wigner surmise. Furthermore, there is no 
indication of spin degeneracy which would result in a bi- 
modal peak spacing distribution |12] ). In Refs. ]lO| , pT[ the 
fluctuations of the peak spacings are considerably larger 



than expected from RMT, whereas the experiments pre- 
sented in Ref. |12j yield smaller peak spacing fluctuations, 
which, however, are still larger than those predicted by 
RMT. 

The deviations from the RMT predictions have been 
frequently interpreted as fluctuations in the charging en- 
ergy. As the charging energy reflects the Coulomb inter- 
actions both between the electrons on the dot as well as 
between the dot and its environment, the dependence of 
the fluctuations on the interaction strength is of funda- 
mental interest. Numerical studies suggest that the fluc- 
tuations are proportional to the charging energy rather 
than to the single particle level spacing |T^,|I^,|l^] . This 
is also found theoretically for the classical limit of a 
Coulomb glass island |l5|]. In RPA calculations the fluc- 
tuations have been related to fluctuations of the eigen- 
functions of the dot jL6) . Another approach based on den- 
sity functional theory emphasizes the role of the Coulomb 
matrix elements of scarred wavefunctions [l7| . Recently, 
also a non-interacting explanation for the Gaussian shape 
of the peak spacing distribution has been given in terms 
of level dynamics due to shape deformation of the quan- 
tum dot @. 

Here we present an experimental study of the statis- 
tics of Coulomb blockade peak positions of a quan- 
tum dot. The dot is defined by patterning the two- 
dimensional electron gas of a silicon metal oxide semicon- 
ductor field effect transistor (MOSFET) structure using 
stacked gates. These experiments differ significantly in 
two major ways from prior experiments on quantum dots 
defined in GaAs/AlGaAs heterostructures: first and fore- 
most, due to the different electron density and material 
properties of silicon, the ratio of the charging energy, Ec, 
to the single particle energy level spacing, Ae, is consid- 
erably larger; likewise the dimensionless parameter r s , 
which characterizes the strength of the Coulomb inter- 
actions, is larger than in previous studies. Secondly, the 
number of electrons is varied by the application of a volt- 
age to a top gate instead of by squeezing the quantum 
dot with a plunger gate. We find that the distribution 
of the peak spacings is unimodal and roughly Gaussian. 
The magnitude of the fluctuations is 15 times larger than 
that predicted by RMT. 

Conduction through a small electron island coupled to 
leads via tunnel barriers is normally suppressed if ksT <C 
Ec, where Eq is the charging energy of the island. This 
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effect is known as the Couiomb biockade m. The block- 
ade is lifted when the condition fid < fidot < Ms is satis- 
fied, where fi s , fid, and fidot are the chemical potentials 
of the source, drain and dot, respectively. In the linear 
response regime, where these experiments have been per- 
formed, \fid — Ms I ^ Ae, IcbT. The chemical potential of 
the dot is defined as fidot(N+l) = E(N+1)-E(N) where 
E(N) is the total energy of the dot occupied by iV elec- 
trons. In the case where the blockade is lifted an electron 
can tunnel from the source onto the dot, changing the 
dot's occupation from N to N + 1, and sequentially tun- 
nel off the dot to the drain leaving the dot in its original 
state. The resulting fluctuation of the electron number 
on the dot leads to a finite conductance. Experimentally 
this can be achieved by appropriately tuning fidot with an 
external gate. Sweeping the gate voltage, V g , results in 
the well known conductance oscillations indicating suc- 
cessive filling of the dot with single electrons. The differ- 
ence of fidot between two adjacent conductance maxima 
is thus given by Afi N = E(N + 1) - 2E(N) + E(N - 1) 
which can be viewed as the discrete second derivative of 
the quantum dot energy with respect to particle number, 
i.e. the inverse compressibility d/i/dN 

In the constant interaction (CI) model 14] the energy 



of the dot is approximated as E(N) = (Ne) /2Cs + 
St=i e «i where the electrostatic interactions are treated 
using a simple capacitive charging model with a total dot 
capacitance Cs, and the quantum mechanical terms are 
taken into account as single particle energies e^. In this 
model the difference of the chemical potentials for succes- 
sive occupation numbers, the so called addition energy, is 
A/ijy = Ec+Aejy with the charging energy Eq — e 2 /C^., 
and the level spacing Aejv = ejv+i — £n- This is mapped 
to gate voltages via e-^AV g = Ec + Aex where C g is the 
capacitance of the dot to the gate and AV g the difference 
between the gate voltages at which adjacent conductance 
maxima occur. 

This final expression motivated the original investiga- 
tions of the peak spacings in the light of random ma- 
trix theory. RMT shows that the normalized spacings S 
((S) = 1) between adjacent eigenvalues of a generic time- 
reversal invariant Hamiltonian are distributed according 
to the Wigner surmise 
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The fluctuations of these spacings are ((S 2 ) — (S*) 2 ) 1 / 2 w 
0.52(5). However, experiments have shown that the 
combined CI-RMT model is not capable of describing 
the observed peak spacing distribution correctly jio] |l2| . 

The Coulomb blockade measurements on which the fol- 
lowing analysis is based have been performed on a quan- 
tum dot defined in a silicon MOSFET structure. We have 
utilized a stacked gate structure to pattern the electron 
gas as shown in Fig. |l|. First, a gate oxide is grown 
on a p-type silicon substrate (lower oxide), and then a 
lower metal gate is deposited and patterned using elec- 
tron beam lithography and lift-off techniques. Above the 



lower gate a second layer of silicon dioxide is deposited 
(upper oxide), and finally an upper gate is formed; the 
upper oxide layer serves to insulate the lower gate from 
the upper gate. Application of positive voltages to the 
upper gate leads to the formation of a two-dimensional 
electron gas (2DEG) at the Si/Si02 interface; n + im- 
planted regions serve as Ohmic contacts to the 2DEG. 
Further details about this device may be found else- 
where 19 2y]. The lower gates locally screen the electric 
field created by the upper gate, and a quantum dot is 
formed by applying appropriate negative voltages to the 
lower gates. The size of the dot is estimated from the 
capacitance to be A as 200 nm x 200 nm, which agrees 
well with the lithographic dimensions of 250 nm x 270 nm 
when electrostatic depletion at the edge is considered. 
The electron density can be varied by changing the up- 
per gate voltage, whereas the lower gate voltage controls 
the tunnel barriers and the electrostatic confinement po- 
tential of the quantum dot. This technique allows the 
definition of very small structures which therefore have 
low capacitances and high charging energies. For the 
quantum dot discussed here these values are Cs ~ 85 aF 
and Ec ~ 1.9 meV as obtained from temperature depen- 
dence measurements of the conductance resonances [ p0| . 

In contrast to previous experiments on quantum dots 
in GaAs/AlGaAs heterostructures the electron density is 
considerably higher, n s 2.5 • 10 16 m~ 2 . The mobility 
of the two-dimensional electron gas is fi = 0.56m 2 /Vs, 
and the mean free path I as 100 nm is comparable to the 
system size. The single particle energy level spacing can 
be obtained from the estimated dot area, A, via Weyl's 
formula Q as Ae = ^t A =15 //eV, where g is the de- 
generacy of electronic states in the two-dimensional elec- 
tron gas, and in* is the effective mass of the electrons. 
In a 2DEG in a silicon MOS system m* = 0.2m e , and 
at B = both the spin and valley degeneracies must be 
considered and therefore g = 4. While these quantum 
dots are smaller than many of the GaAs/AlGaAs quan- 
tum dots studied |]l0|-|l2]| , Ae is of the same order due to 
the larger effective mass, and to the valley degeneracy. 

The strength of the electronic interactions char- 
acterized by the dimensionless parameter r s = 
gj (2y/Trn s a* B ) — 2.1 considerably exceeds the values ob- 
tained in recent experiments (where r s ral) 10-^]; here 
a* B is the effective Bohr radius. Similarly, the ratio of 
charging energy to single particle energy level spacing 
Eq / Ae f=a 125, another measure of the relative impor- 
tance of electron-electron interactions, is larger than in 
the experiments performed on GaAs/AlGaAs quantum 
dots. 

The measurements were performed in a 3 He refriger- 
ator at a temperature of T = 320mK using standard 
lock-in techniques at low frequencies and bias. The con- 
ductance oscillations were measured as a function of the 
upper gate voltage. Consequently, the electron density 
was varied without drastically changing other system pa- 
rameters such as charging energy, single particle energy, 
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and dot shape. This also contrasts with former experi- 
ments on the statistics of conductance oscillations where 
the shape of the quantum dot was distorted by plunger 
gates gH- 

The following analysis is based on a series of more 
than 100 conductance peaks occurring in the upper gate 
voltage range from 12.1 V to 13.5 V (see Fig. || in- 
set). In this range the electron density changes from 
2.4 x 10 16 m~ 2 to 2.6 x 10 16 m~ 2 . The quantum point 
contacts connecting the quantum dot to the leads are 
tuned into the tunneling regime by applying voltages of 
—4.5V and —8.0V to the left and right pair of lower 
gates, respectively. The position of each peak is obtained 
by fitting the peak by a thermally broadened line shape 
oc cosh~ 2 (ea(V — Vb)/2fc_BT) j| where a and Vb serve as 
the fit parameters. The gate voltage spacings AVg are 
calculated from the peak positions. 

The mean value (AVg) of the spacings as a func- 
tion of V g is obtained from a linear fit as AV g (V g ) sa 
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x (T/ ff - 12.5 V) 



mV. The smallness of 
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the slope of this fit shows that the influence of the upper 
gate on the capacitance and therefore on the size of the 
dot is rather weak. Accordingly, the shape deformation 
which has been postulated to explain the distribution of 
AVg p8[ | plays no significant role in this experiment. The 
normalized peak spacings 
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are displayed in Fig. |2|. The fluctuations of 8 
puted to be (8 2 ) 1/2 « 0.06. The fluctuat ions in the ad- 
dition energy are, therefore, roughly 115/ieV, which is 
7.5 times the mean level spacing Ae, and thus 15 times 
larger than the fluctuations expected from RMT. This 
supports the view that the fluctuations of the addition 
energy scale with the Coulomb energy rather than with 
the kinetic energy. However, the proportionality factor 
0.06 is smaller than that suggested by numerical calcu- 
lations (0.1 - 0.2) |Q. It should be noted that in these 
experiments Ae w fc^T. We expect that the effect of 
thermal broadening would be to reduce the fluctuations 
in peak spacing. A simple model |l2j predicts that the 
fluctuations expected within RMT would be reduced by 
a factor of 2 - 3. If we incorporate this correction into 
RMT then the fluctuations we find in our experiment are 
30 - 45 times larger than those predicted by RMT. 

The distribution of the peak spacings normalized to 
an area of unity is shown in Fig. |3|. The distribution is 
unimodal and roughly has the shape of a Gaussian. In the 
inset of Fig. ^ the experimental distribution is depicted 
together with the Wigner surmise (Eq. [l]) ; for comparison 
to 8 we have rescaled the predictions of RMT taking into 
account the experimental values of Eq and Ae. As in 
previous experiments there is no evidence of a bimodal 
addition spectrum as is predicted by the CI- RMT model. 
This is in agreement with the theoretical prediction that 



the influence of spin degeneracy on the addition spectrum 
is washed out for stronger electron-electron interactions 
(r. > 1) @0. 

In conclusion, we have investigated the Coulomb block- 
ade peak spacing distribution of a quantum dot fabri- 
cated in the 2DEG of a silicon MOSFET structure. In 
accordance with experiments on GaAs/AlGaAs quantum 
dots the distribution differs from the Wigner surmise and 
is roughly Gaussian. The fluctuations are approximately 
0.06 x Ec- Due to the large ratio of charging energy Eq 
to single particle energy Ae this strongly suggests that 
the fluctuations scale with Eq and not with Ae. This 
clearly contradicts the predictions of RMT and indicates 
that the fluctuations are dominated by electron-electron 
interactions in this system. 
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FIG. 1. Schematic representation of the device design. A 
cross section of the device is shown in (a). Two oxide and two 
gate layers are formed on top of p-type silicon substrate. The 
lower oxide layer has a thickness of roughly 20 nm, while the 
upper oxide is approximately 80 nm thick. The voltage on the 
upper gate is used to vary the electron density in the 2DEG 
induced at the interface of the lower oxide and the silicon. A 
top view of the device is shown in (b). The pattern in the 
lower gates defines a quantum dot in the induced electrons; 
note that the upper gate covers all of the area show in (b) and 
overlaps the source and drain. The lithographic dimensions 
of the quantum dot are 250 nm x 270 nm. 

FIG. 2. The normalized peak spacings 8 obtained from an 
upper gate voltage sweep. The fluctuations around the mean 
value are much larger than expected from the CI-RMT 
model. The inset shows Coulomb blockade conductance os- 
cillations as a function of the upper gate voltage. 



FIG. 3. Histogram showing the distribution of the normal- 
ized peak spacings from Fig. 2. The area of the histogram 
is normalized to unity. A Gaussian fit with standard devi- 
ation of a = 0.06 is also displayed. The inset shows the 
same histogram alongside the Wigner surmise, the distribu- 
tion predicted by RMT, for Ae = 15 fj,eV. The experimental 
distribution is much broader than expected from RMT. 
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